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Abstract. We summarize the higher-loop perturbative computation of the ghost and gluon propagators in SU (3) Lattice 
Gauge Theory. Our final aim is to compare with results from lattice simulations in order to expose the genuinely non- 
perturbative content of the latter. By means of Numerical Stochastic Perturbation Theory we compute the ghost and gluon 
propagators in Landau gauge up to three and four loops. We present results in the infinite volume and a — > limits, based on 
a general fitting strategy. 
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1. INTRODUCTION 

This talk summarizes our work on the higher-loop pertur- 
bative gluon and ghost propagators in Landau gauge [1, 
2]. The Monte Carlo study of both propagators, which 
are closely related to each other by Schwinger-Dyson 
equations (SDE), has attracted much attention outside the 
lattice community by phenomenologists working on in- 
frared QCD in general and hadron physics (see our orig- 
inal papers for further references). Taken together, both 
propagators provide us with a definition and the momen- 
tum dependence of the running coupling (Xs{q^) directly 
based on the ghost-gluon vertex. 

A simple connection between the two propagators 
exists in the extreme infrared, both beeing powerlike 
in a scaling or massive in a decoupling solution. This 
nonuniqueness reflects the Gribov problem. The effect 
of nontrivial vacuum structure (vortices, instantons) is 
manifest also in the gluon propagator, in the intermedi- 
ate momentum range around ^(1 GeV) where the SDE 
approach suffers from truncation ambiguities and where 
nonperturbative lattice calculations are unrivalled. In or- 
der to follow the onset of nonperturbative effects, it is 
desirable to approach this momentum range from high 
momenta within higher-order perturbation theory. While 
ordinary diagrammatic lattice perturbation theory (LPT) 
soon gets too involved to be pursued. Numerical Stochas- 
tic Perturbation Theory (NSPT, for a recent review see 
Ref. [3] and references therein), provides a powerful tool 
to perform high-loop computations. 
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2. NSPT IN A NUTSHELL 

NSPT has its roots in stochastic quantization and is based 
on a modified Langevin equation equipped with stochas- 
tic gauge fixing. We use here a version for quenched lat- 
tice QCD with Wilson gauge action. Actually, it is a hi- 
erarchy of first-order evolution equations associated with 
various parts of the gauge link fields U and gluon fields 
A exposed by an expansion in powers of the lattice cou- 
phnggoc 1/y^: 

f/.. = irM^,A,i,, = lr'/^A« (1) 

These different orders are separately dealt within the 
code. The maximal addressable order of perturbation the- 
ory is thus limited by the available computing resources 
(cpu time and memory). 

The Langevin simulation is implemented in an Euler 
scheme with a finite evolution time step. Before the es- 
timator for the gauge dependent ghost and gluon propa- 
gators can be evaluated, we have to fix the gauge to the 
minimal Landau gauge. For this purpose, a sequence of 
configurations (separated by i^(50) Langevin time steps) 
is subjected to a Fourier-accelerated gauge-fixing proce- 
dure, after which the individual gluon fields, aJ^'^ (associ- 
ated with particular perturbative order g') are transversal 
within machine precision. 

The propagators are evaluated taking the long-time av- 
erage of coefficients, order by order in a loop expansion 
in even powers of g. Contributions from odd powers van- 
ish within the statistical errors. As for any Langevin sim- 
ulation, one then has to take the limit to vanishing time 
step. In order to get results comparable with the prac- 



tice of LPT, the continuum limit and the limit of infinite 
volume must be performed. NSPT results for finite lattice 
volume and spacing can be confronted directly with stan- 
dard MC results for a given j3, provided the definitions of 
the studied observables is the same. 

The gluon two-point function in n-loop order is de- 
fined as a convolution of the bilinears of gluon fields 
(in momentum space) in complementary orders (A^l' = 
T"A<'\p^{k^)^27:k^/{aN)y. 
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In Landau gauge we consider = 3Z)'"^ and use 

the dressing functions (pn{k^) = {2/ a) sin [izk^/N)) 

J^^^"\p) = / D''"\p{k)),P-(''\p) = p2 D''"\p{k)) . 

(3) 

The color diagonal ghost propagator in momentum 
space is the color trace in the adjoint representation 



G{p{k)) = -{Tu^,M-\k))^ 



(4) 



In (4) M^^{k) is the Fourier transform of the inverse 
FP operator in lattice coordinate space. It is expanded 
in terms of products of various A^'\ with the term M'-'^ 
collecting all terms of order . This structure allows to 
express M^' [k) also as an expansion in orders of g in a 
recursive way. Again we use the ghost dressing functions 

=p2G(")(/,(fc)), =/32dW(p(^)). (5) 

As an example the cumulatively summed perturbative 
gluon dressing function for various volumes is shown in 
Fig. 1 using — 6/15 ~l. 
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FIGURE 1. The cumulatively summed perturbative gluon 
dressing function for various volumes. 

A reasonable "convergence" of the NSPT results up 
to few loops (three or four are available now) requires a 
small bare coupling g. However, g is known to be a poor 
expansion parameter [4]. One can speed up convergence 



by "boosting", i.e., trading the bare coupling constant by 
an effective "boosted" coupling g| = / Ppenig^) > g^- 
Here Ppert is defined by the average perturbative plaque- 
tte determined also within our Langevin simulations. The 
effect of the boosted coupling being larger is overcom- 
pensated by the rapid decay of the expansion coefficients 
with increasing order n. 

We illustrate the effect of "boosting" the perturbative 
expansion and confront the boosted dressing functions 
with corresponding new Monte Carlo (MC) data of the 
Berlin group [5] adopting the same definitions for the 
propagators and the gauge fixing as in NSPT. This is 
shown in Figs. 2,3 where also the bare and the boosted 
inverse couplings /3 and jSboost given. As expected, 
boosting moves the NSPT data closer to the MC results, 
but they cannot be reached completely, certainly not at 
j3 = 6.0. 
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FIGURE 2. The ghost dressing function in three loops for 
naive and boosted NSPT comapred to MC data for a 12'* lattice. 
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FIGURE 3. Same as in Fig. 2 for the gluon dressing function 
in four loops. 

Here we define the renormalization-group invariant 
running coupling by the ghost-gluon vertex in a par- 
ticular (minimal) MOM scheme (see e.g. [6]). It is given 
in terms of the bare gluon and ghost dressing functions 
and J as follows: 



J{p{k),l5fj^{p{k),P) 



4;rj3 



(6) 



The a.s calculated from the NSPT dressing functions, 
both summed up to the orders available, is compared to 
the MC results at )3 = 6.0 and 9.0. The corresponding 



data is shown in Figs. 4,5 again for naive and boosted 
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FIGURE 4. Comparing naive and boosted PT (based on 
NSPT data) for the running coupling constant to correspond- 
ing MC data for a 12'* lattice at jS = 6.0. 
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FIGURE 5. Same as in Fig. 4 at j6 = 6.0 zoomed into the 
large momentum region. 

perturbation theory. We see that the running coupling 
from MC simulations is approached at large momenta 
from below up to 7% for j3 = 6.0 and practically ap- 
proached within the present errors for j3 = 9.0, i.e. in an 
effectively deconfined phase. 

3. DRESSING FUNCTIONS IN THE 
LIMITS V^oc AND ap^O 

In the RI'-MOM scheme the dressing functions are just 
the wave function renormalization constants for the ghost 
and gluon fields at the renormalization point jx^ — p^. 
The NSPT data available at various volumes and lattice 
momentum realizations allow us to find the perturbative 
dressing functions to three-loop accuracy in the bare cou- 
pling including the non-logarithmic contributions. Via 
standard transformations the results can be transformed 
to the renormalized coupling in the preferred scheme. To 
find those constants, a fitting procedure has been pro- 
posed which takes into account both hypercubic and fi- 
nite volume effects (for details see [1]). As result we get 
(L = log(pfl)2) 



J{a,p,l5) = 1 



-0.0854897L + 0.525314 
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-8.13(3)L + 31.7(5) 



0.08211l2-1.48445L + 7.93(12) 



-0.02964L' +0.81689 



(8) 



The results for one-loop lattice perturbation theory are 
known for a long time, the higher-loop non-leading log's 
and constant contributions are our predictions for the 
Landau gauge. 

4. CONCLUSION 

We have calculated the gluon propagator in Landau 
gauge up to four and the ghost propagator up to three 
loops in NSPT. The dressing functions summed using 
boosted PT are compared to recent MC measurements 
of the Berlin Humboldt University group. At large lat- 
tice momenta the dressing functions with more than four 
loops will match the MC measurements, thus enabling 
a fair accounting of the perturbative tail taking care of 
discretization effects. This can be used as an alterna- 
tive to fitting the high momentum tail of MC results by 
continuum-like formulae. 

We worked out the relation to standard LPT in limits 
y — > oo and pa — > 0. For this aim we developed a fitting 
strategy for lattice artifacts and finite-size corrections. 
We find good agreement with known one-loop results of 
diagrammatic LPT and present original two- and three- 
loop results for the propagators. 
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